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Attached  is  a  Fortran  subroutine  (TIMESQ)  and  a  companion  subroutine 
(GADSSPN)  to  generate  essentially  a  continuous  clutter  signal  sequence  where 
the  parameters  are  set  in  real-time.  The~mgtttPd~±s~baaed  on  that  in  Reference  ' 
"1.  For  this  preliminary  version,  three  arrays  specifying  the  amplitude 
spectral  functions  on  the  JI^AA,  and  AE  channels  are  the  inputs  to  TIMESQ. 

^  _  The subroutine  to  generate  the  amplitude  spectral  functions  from 

the  engagement  geometry  is  under  development.  Basically  a  table  of  parameter 
values  will  be  stored;  then  given  the  engagement  geometry  the  desired  values 
will  be  retrieved  from  the  table,  possibly  with  interpolation;  finally,  the 
amplitude  spectral  functions  will  be  generated  from  the  parameters.  Current 
design  is  based  on  14  parameteters  for  the  three  channels.  The  most  time 
consuming  computation  predicted  is  an  exponential  (4  ysec  on  the  AP120B) 
for  each  spectral  sample.  The  portion  of  the  FFT  computation  for  each 
sample  is  6.3  ysec,  and  the  total  time  allowed  for  all  operations  is  27 
ysec  per  sample.  Therefore,  there  appears  to  be  sufficient  time  available 
to  implement/ this  approach  on  the  AP120B. 
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SUBROUTINE  TIMESQ (SA> 

C 

C  IN  THIS  SUBROUTINE  WE  GENERATE  A  CORRELATED  TIME  SEQUENCE  THAT  HAS  AN 
C  AMPLITUDE  SPECTRAL  FUNCTION  DESCRIBED  BY  ARRAY  SA  OF  LENGTH  NFFT.  THE 
C  METHOD  USED  IS  THE  ON-LINE  FFT  APPROACH  DESCRIBED  IN  SECTION  6  OF  MR I 
C  REPORT  132-44,  WHERE  SUCCESSIVE  CALLS  TO  TIMESQ  WILL  CREATE  CONTIGUOUS 
C  TIME  SEQUENCES  BASED  ON  OVERLAPPING  FFTS  IN  THE  FREQUENCY  DOMAIN. 

C 

C  ON  INPUT . 

C 

C  NFFT  =  SIZE  OF  FFT  < =2** INTEGER ) 

C  NOVLP  =  NUMBER  OF  SAMPLES  THAT  THE  FFTS  OVERLAP  (NFFT/B  IS  OK) 

C  SA(K)  -  SPECTRAL  AMPLITUDE  OF  KTH  SAMPLE,  K»1,...,NFFT 

C_ _ __  _ _ 

C  ON  OUTPUT . 

C 

C_ _ XR<K)  -  KTH  SAMPLE  OF  IN-PHASE  SIGNAL,  K«1,...,NF 

C  XI(K)  -  KTH  SAMPLE  OF  QUADRATURE  SIGNAL,  K»1,...,NF 

C 

C  WHERE  NF-NFFT -NOVLP . _ 

C 

C  THE  ARRAYS  MUST  BE  DIMENSIONED  AS  LARGE  AS . 

C _ 

C  SA  NFFT 

C  XR  NFFT +NOVLP— 1 

C _ XI  NFFT+NOVLP-1 

C  A  NOVLP -1 

C 

C  THIS  SUBROUTINE  HAS  NOT  BEEN  TESTED 
C 

DIMENSION  SA ( 1 ) , A ( 64 ) 

COMMON  /C2/_  NFFT , NOVLP 
COMMON  /C4/  XR  <  320 ) , X I ( 320 ) 

DATA  NO/O/ 

IF  (NO.  GT.  O)  GO  TO  15 
NO-1 

NF-NFFT-NOVLP 
Nl— NOVLP— 1 
DO  10  K— 1, Nl 
A ( K ) -K/FLOAT  <  NOVLP ) 

XR  <  K+NFFT )  -0. 

XI <K+NFFT)-0. 

10  CONTINUE 

15  CALL  OAUSSPN(SA,  XR,  XI,  NFFT) 

CALL  FFT2(XR, XI, NFFT, 1 ) 

DO  20  K— 1,  Nl 

_  XR  < K ) -A ( K ) *  <  XR ( K ) -XR  <  K+NFFT ) ) +XR <  K+NFFT ) 

XI(K>-A(K)*(XI(K)-XI (K+NFFT) )+XI (K+NFFT) 

XR (K+NFFT ) -XR ( K+NFFT-NOVLP ) 

_XX ( K+NFFT ) -X I ( K+NFFT-NOVLP ) 

„  30  CONTINUE 
. RETURN 
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SUBROUTINE  CAUSSPN-(A,  XR,  XI#  N) 

C  .  .  . 

C  OENERATES  RANDOM  PHASOR  COMPONENTS  < XR(K)>  XI (K) )  OF  ZERO  MEAN  WITH  THE 
C  AVERAGE  POWER  OF  XR<K)**2+XI <K>**2  GIVEN  BY  A<K)**2  FOR  K=1,...,N. 

C  _ 

C  THE  RANDOM  NUMBERS  ARE  ACCESSED  FROM  RECIRCULATING  TABLES.  NTR  AND 
C  NTI  ARE  THE  CIRCULATION  PARAMETERS  AND  NMAX=MAX  <N) .  WE  SHOULD  CHOOSE 
C  NTR  AND  NTI  RELATIVELY  PRIME. 

C 

C  ARRAYS  TR  AND  TI  MUST  BE  DIMENSIONED  AS  LAROE  AS  NTR+NMAX  AND  NTI+NMAX 

C  RESPECTIVELY.  •_  _ 

C 

C  GAUSS  IS  A  FUNCTION  SUBROUTINE  THAT  GENERATES  A  GAUSSIAN  R.  V.  OF  ZERO 

C  MEAN  AND  UNIT  VARIANCE. _ _ _  _ _ 

C  MiT 

C  IF  IRSET-1  IN  COMMON  /GS/>  THE  POINTERS  WILL  BE  RESTORED *T0  THEIR 
__£_J3RI0INAL_VALUES  WHEN  GAUSSPN  WAS  KMTtftfcP.  THIS  MEANS  THAT  THE  SAME 
C  GAUSSIAN  SEQUENCE  (EXCEPT  FOR  THE  WEIGHTS  A>  WILL  BE  GENERATED  ON  THE 
C  NEXT  CALL. 

_ : _ . -  _ 

C  THIS  SUBROUTINE  HAS  NOT  BEEN  TESTED 
C 

_ DIMENSION  A<1>. XR(1>. XI<1) 

DIMENSION  TR< 1253).  TI ( 1353) 

COMMON'  /GS/  IRSET 

DATA  NTR/ NTI/997*  1097A  NMAX/256/,  KRi  KI/O*  1/ 

IF(KR.  OT.  0)  GO  TO  30 
IF(N.  OT.  NMAX)  STOP 

_ Nl-NTR+NMAX 

DO  20  k*l.Nl 

TR(K)«.  707 107*CAUSS( DUMMY) 

_  20. CONTINUE 

N2-NTI+NMAX 
DO  25  K*l<  N2 

_  TICK)-.  707107*GAUSS< DUMMY) 

25  CONTINUE 
KR-1 

30  DO  35  I-i.N  _ 

XR( I )— A( I >*TR(KR ) 

XI ( I )— A< I)*TI(KI) 

KR-KR+1 
KI-4<i+l 
35  CONTINUE 

_IFC  IRSET.  NE.  1)  00  TO  40 
KR-KR-N 
KI-KI-N 

40  JFCKR.  OT.  NTR)  KR-KR-NTR 
IFCKI.GT.  NTI)  KI-KI-NTI 
RETURN 

_ END _ 
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